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Calculation of proton-^He elastic scattering between 7 and 35 MeV 
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Background: Theoretical calculations of the four-particle scattering above the four-cluster breakup threshold 
are technically very difficult due to nontrivial singularities or boundary conditions. Further complications arise 
when the long-range Coulomb force is present. 

Purpose: We aim at calculating proton-^He elastic scattering observables above three- and four-cluster breakup 
threshold. 

Methods: We employ Alt, Grassberger, and Sandhas (AGS) equations for the four-nucleon transition operators 
and solve them in the momentum-space framework using the complex-energy method whose accuracy and practical 
applicability is improved by a special integration method. 

Results: Using realistic nuclear interaction models we obtain fully converged results for the proton-'^He elastic 
scattering. The differential cross section, proton and ^He analyzing powers, spin correlation and spin transfer 
coefficients are calculated at proton energies ranging from 7 to 35 MeV. Effective three- and four-nucleon forces 
are included via the explicit excitation of a nucleon to a A isobar. 

Conclusions: Realistic proton-'^He scattering calculations above the four-nucleon breakup threshold are feasible. 
There is quite good agreement between the theoretical predictions and experimental data for the proton-^He 
scattering in the considered energy regime. The most remarkable disagreements are the peak of the proton 
analyzing power at lower energies and the minimum of the differential cross section at higher energies. Inclusion 
of the A isobar reduces the latter discrepancy. 

PACS numbers: 21.45.-v, 25.10.-l-s, 21.30.-x, 24.70.-l-s 
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I. INTRODUCTION 

Froton-'^He (p-'^He) scattering is one of the most com- 
monly used experiments to study the four-nucleon sys- 
tem [ij: It involves two charged particles that are sta- 
ble and easy to detect with acceptable precision, there 
are no competing channels until Ep — 7.3 MeV proton 
lab energy, and beyond that only three- and four-cluster 
breakup takes place up to the pion-production thresh- 
old. Much like neutron-"^H (n-'^H) scattering, p-'^He is 
dominated by isospin T = 1, but has the deciding ex- 
perimental advantage of having a proton beam and a 
non-radioactive '^He target. On the contrary, from the 
theoretical point of view, p-^He is more difficult to cal- 
culate than n-'^H due to the long-range Coulomb force 
between protons (p) that gives rise to complicated bound- 
ary conditions in the coordinate-space and non-compact 
kernel in the momentum-space representation. Neverthe- 
less, these difficulties have been solved below the three- 
cluster breakup threshold using three different theoretical 
frameworks, namely, the hyperspherical harmonics (HH) 
expansion method 0, Q , the Faddeev-Yakubovsky (FY) 
equations [3 for the wave function components in coordi- 
nate space [a, Q, and the Alt, Grassberger and Sandhas 
(AGS) equations Q for transition operators that were 
solved in the momentum space [1, Q- A good agree- 
ment between these methods has been demonstrated in a 
benchmark for n-^H and p-^He elastic scattering observ- 
ables [lOJ using realistic nucleon-nucleon (NN) potentials. 

Recently we extended the AGS calculations to energies 
above three- and four-cluster breakup thresholds pTlll3|. 
The complex energy method [IJ, [l4|] was used to deal 



with the complicated singularities in the four-particle 
scattering equations; its accuracy and practical appli- 
cability was greatly improved by a special integration 
method [il|. This allowed us to achieve fully converged 
results for n-'^H elastic scattering and neutron-neutron- 



deuteron recombination into 



^H using realistic NN 



interactions. We note that the FY calculations of n-H 
elastic scattering have been recently extended as well to 
energies above the four-nucleon breakup threshold [15j . 
however, using a semi-realistic NN potential limited to 
S"- waves. 

In the present work we extend the method of Ref. [Il| 
to calculate the p-'^He elastic scattering above breakup 
threshold and compare with existing data for cross sec- 
tions and spin observables over a wide range of proton 
beam energies up to Ep — 35 MeV. The pp Coulomb in- 
teraction is treated as in Refs. 0, Ha] using the method 
of screening of the pp Coulomb potential followed by the 
phase renormalization of transition amplitudes [ITI . llSj . 
Thus, standard AGS scattering equations with short- 
range potentials are applicable. At energies bellow three- 
cluster threshold our results agree with those obtained 
by other methods as mentioned in Ref. 10]. Compared 
to our previous n-'^H scattering calculations above the 
breakup threshold [ll| , the most serious complication for 
p-'^He is the convergence of the partial-wave expansion 
that requires a larger number of states due to the longer 
range of the screened Coulomb potential. 

In Sec.[TT]we describe the theoretical formalism and in 
Sec, mil we present the numerical results. The summary 
is given in Sec. |TVl 



II. 4N SCATTERING EQUATIONS 

We use the symmetrized AGS equations @ as appro- 
priate for the four-nucleon system in the isospin formal- 
ism. They are integral equations for the four-particle 
transition operators Upa, i-e., 

+ U2GotGoU2u (la) 

i^2i = (Go t Go)-\l - F34) + (1 - PsiWiGo t GoUn- 

(lb) 

Here, a = 1 corresponds to the 3 -f 1 partition (12,3)4 
whereas a = 2 corresponds to the 2-|-2 partition (12) (34); 
there are no other distinct two-cluster partitions in the 
system of four identical particles. The free resolvent at 
the complex energy E + is is given by 



Go^iE + ie^ Ho) 



(2) 



with Hq being the free Hamiltonian. The pair (12) tran- 
sition matrix 



t — V + vGot 



(3) 



is derived from the potential v; for the pp pair v includes 
both the nuclear and the screened Coulomb potential wr. 
Our calculations are done in momentum space; however, 
we start with the configuration-space representation 



WR{r)^wir)e~'-''^'^^^' 



(4) 



where w{r) = a^/r is the true Coulomb potential, a^ ~ 
1/137 is the fine structure constant, R is the screening 
radius, and n controls the smoothness of the screening. 
All transition operators acquire parametric dependence 
on R but it is suppressed in our notation, except for 
the scattering amplitudes. The symmetrized 3-1-1 or 2-1-2 
subsystem transition operators are obtained from the re- 
spective integral equations 



Ua — Ja^n 



Pat Go Ua 



(5) 



The basis states are antisymmetric under exchange of two 
particles in the subsystem (12) for the 3 + 1 partition 
and in (12) and (34) for the 2 -I- 2 partition. The full 
antisymmetry of the four-nucleon system is ensured by 
the permutation operators Pat of particles a and b with 

Pi = P12 P23 + Pi3 P23 and P2 = Pi3 P2A- 

The p-^He scattering amplitude with nuclear plus 
screened Coulomb interactions at available energy E = 
ei -|- 2pi/3mpf is obtained from the on-shell matrix ele- 
ment (pi|T(fl.)|pi) — 3{(l)'i\Uii\(f>i) in the hmit e — > +0. 
Here |0i) is the Faddeev component of the asymptotic 
p-^He state in the channel a = 1, characterized by the 
bound state energy ei = —7.72 MeV and the relative 
p-^He momentum pi, niN being the average nucleon 
mass. Due to energy conservation p[ = pi. 

The amplitude (p']^|T(^-)|pi) is decomposed into its 
long-range part (p']^|i^™|pi), being the two-body on-shell 



transition matrix derived from the screened Coulomb po- 
tential of the form ^ between the proton and the cen- 
ter of mass (cm.) of "^He, and the remaining Coulomb- 
distorted short-range part. Renormalizing (p']^|T(fl)[pi) 
by the phase factor Z^^ [^, [13, [3 1 i'^ the i? — ?► 00 limit, 
yields the full p-^He transition amplitude in the presence 
of Coulomb 



(pl|T|pi) = (pl|i^-|pi) 



+ hm {(p;|[r(^)-ir]|Pi)^fi'}, 



(6) 



where the first term is obtained from Z^^(p']^|t^™|pi) 
that converges, in general as a distribution, to the ex- 
act Coulomb amplitude (pi |t"°ipi) between the proton 
and the cm. of the ^He nucleus [13. The renormaliza- 
tion factor Zr is defined in Refs. 9, 16]. The second 
term in Eq. ([6]) , after renormalization by Zjj , represents 
the Coulomb-modified nuclear short-range amplitude. It 
has to be calculated numerically, but, due to its short- 
range nature, the i? — > cxd limit is reached with sufficient 
accuracy at finite screening radii R. Since the conver- 
gence with R is faster at higher energies, the required 
screening radii are smaller than in our low-energy p-^He 
calculations Q. We found that R ranging from 8 to 10 
fm leads to well-converged results in the energy regime 
considered in the present paper. Furthermore, we take 
a sharper screening with n = 6 such that at short dis- 
tances r < R the screened Coulomb approximates the 
full Coulomb better than with n = 4 used in Ref. Q and 
at the same time vanishes more rapidly at r > i? thereby 
accelerating the partial- wave convergence. 

We solve the AGS equations ([T]) in the momentum- 
space partial-wave framework. The states of the total 
angular momentum J with the projection A4 are 
defined as \kxkykz[lz{{ly[{lxSx)jxSy\Sy}JyS^)Sz\JM) 
for the 3-1-1 configuration and 

\kxkykz{lz{{lxSx)jx[lv{sySz)Sy\jy}Sz)JM) for the 
2 -f- 2. Here k^, ky and k^ are the four-particle Jacobi 
momenta in the convention of Ref. |l9| . Ix, ly, and Iz 
are the associated orbital angular momenta, jx and 
jy are the total angular momenta of pairs (12) and 
(34), Jy is the total angular momentum of the (123) 
subsystem, Sy and Sz are the spins of nucleons 3 and 4, 
and Sx, Sy, and Sz are channel spins of two-, three-, 
and four-particle system. A similar coupling scheme 
is used for the isospin. We include a large number of 
four-nucleon partial waves, up to lx,lyjz,jx,jy — 7, 
Jy — ^, and i7 = 6, such that the results are well 
converged. In fact, lower cutoffs are sufficient for lower 
J", e.g., lxjyjz,jx,jy £ 5 and Jy < | are sufficient for 
J' < 3. Furthermore, for most observables J' < 5 or 
even J < 4 are enough for the convergence; J^ — 6 yields 
small but still visible effect only above Ep — 30 MeV. 

The numerical calculations are performed for complex 
energies, i.e., with finite e. The limit e — >■ -1-0 needed 
for the calculation of the amplitude (pi|T(/j)|pi) is ob- 
tained by the extrapolation of finite e results as proposed 
in Ref. 13[. A special integration method developed in 
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FIG. 1. (Color online) Differential cross section for elastic 
p-^'He scattering at 8.52, 13.6, 19.4, 25.0, and 35.0 MeV pro- 
ton energy as function of the cm. scattering angle. Results 
obtained with INOY04 (solid curves), and, at selected ener- 
gies, with CD Bonn (dashed-dotted curves) and AV18 (dotted 
curves) potentials are compared with the experimental data 
fromRefs. ^i^. 



Ref. [Ul is used to treat the quasi-singularities of the 
AGS equations (HJ. We obtain accurate results by using 
£ ranging from 1 to 2 MeV at lowest considered energies 
and from 2 to 4 MeV at highest considered energies. Grid 
points for the discretization of each momentum variable 
range from 30 (at lower energies) to 35 (at higher ener- 
gies). Further details on the other numerical techniques 
for solving the four-nucleon AGS equations can be found 
inRef. [fi]. 



III. RESULTS 

We study the p-^He scattering using realistic high- 
precision NN potentials, namely, the Argonne (AV18) 
potential [20|, the inside-nonlocal outside- Yukawa 
(INOY04) potential by Dolcschall jl, [2l[, the charge- 
dependent Bonn potential (CD Bonn) [22], and its ex- 
tension CD Bonn + A 23] allowing for an excitation of a 
nucleon to a A isobar and thereby yielding effective three- 
and four-nucleon forces (3NF and 4NF). The '^He bind- 
ing energy calculated with AV18, CD Bonn, CD Bonn + 
A, and INOY04 potentials is 6.92, 7.26, 7.54, and 7.73 
MeV, respectively; the experimental value is 7.72 MeV. 
Therefore most of our predictions correspond to INOY04 
as it is the only potential that nearly reproduces the ex- 
perimental binding energy of '^He. The calculations with 
other potentials are done at fewer selected energies. 

In Figs. [T] - [5] we show the differential cross section 



FIG. 2. (Color online) Same as in Fig.[T]but at 10.77, 16.23, 
21.3, and 30.0 Me V p roton energy. The experimental data 
are from Refs. [HM. 



da/dil for elastic p-^He scattering as a function of the 
cm. scattering angle 0c.m. at a number of proton en- 
ergies ranging from Ep = 8.5 to 35.0 MeV. This ob- 
servable decreases rapidly with the increasing energy and 
also changes the shape; the calculations describe the en- 
ergy and angular dependence of the experimental data 
fairly well. Below £^p = 15 MeV the experimental data 
are slightly underpredicted at forward angles as happens 
also at energies below the three-cluster breakup thresh- 
old [3, [i3- At the minimum the da/dfl predictions scale 
with the ^He binding energy: the weaker the ^He bind- 
ing the lower the dip of da/dO. that is located between 
©cm. = 105° and 6c.m. = 125°. The scahng is more pro- 
nounced at higher Ep. For the INOY04 potential that fits 
the '^He binding energy, one gets an good agreement in 
the whole angular region up to Ep ~ 20 MeV but, as the 
energy increases, the calculated cross section underpre- 
dicts the data at the minimum much like what happens 
in nucleon-deuteron elastic scattering [1^, 123, 1231 but for 
nucleon energies above 60 MeV. In line with the conjec- 
tures that were made 15 years ago for the three- nucleon 
system l2a l2Ml ' this underprediction of the data at the 
minimum of da/dQ may be a sign for the need to include 
the 3NF. 

In Fig. [3] we show the proton analyzing power Ay for 
elastic p-'^He scattering at proton energies ranging from 
7.0 to 35.0 MeV. We observe that the sensitivity to the 
nuclear force model and energy is considerably weaker 
as compared to the regime below three-cluster threshold 
13, 13 . Most remarkably, in contrast to low energies where 
the famous p-'^He ^^-puzzle exists @, HH], the peak of 
Ay around 120 degrees is described fairly well but there 
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FIG. 3. (Color online) Proton analyzing power Ay for elastic p-^He scattering at 7.03, 8.52^ 10.03, 13.6, 19.4, 21.3, 25.0, 30.0, 
and 35.0 MeV proton energy. Curves as in Fig. [l] The experimental data are from Refs. [3CH3a |. 



is a discrepancy in the niinimum region. This is similar 
to the energy evolution of the Ay-puzzle in the three- 
nucleon system [35] . 

In Fig. |4] we show the '^He analyzing power Aoy for 
elastic p-^He scattering at proton energies ranging from 
7.0 to 35.0 MeV. Aoy varies slowly with energy but is 
slightly more sensitive to the NN potential. Contrary to 
Ay , calculated AQy is in better agreement with data over 
the whole energy range, particularly when the INOY04 
potential is used. 

The experimental data are scarcer for double polar- 
ization observables. In Fig. [5] we show the p-^He spin 
correlation coefficient Ayy for elastic p-^He scattering at 
7.03, 8.52, 10.03, and 19.4 MeV, and in Fig. [U we show 
Axx for elastic p-^He scattering at 8.52 and 19.4 MeV 
proton energy. Calculated Ayy exhibits some sensitivity 
to the NN potential model and describes the data rea- 
sonably well; the agreement with data is the best when 
the INOY04 interaction is used. The same happens for 
Axx but for the single data set we know of. 

Finally in Fig. [7] we show the proton spin transfer co- 
efficients K^ , K^ , and K^ for elastic p-'^He scatter- 
ing at 8.52, 10.77, and 16.23 MeV proton energy. Note 



that 8.52-MeV predictions are compared to experimen- 
tal data taken at 8.82 MeV but, given the weak energy 
dependence of these observables, the comparison is ap- 
propriate. The calculated spin transfer coefficients show 
a rich angular structure and follow the data reasonably 
well but cannot be fully tested by the available data con- 
fined to the angular region below 0c.m. — 110°. In con- 
trast to other shown spin observables, the spin transfer 
coefficient K'^ around 8c.m. — 120°, i.e., in the region 
of the differential cross section minimum, varies quite 
rapidly with the energy. There is little sensitivity to NN 
interaction model, except for K^ around 8c.m. = 120° at 
Ep — 16.23 MeV. Theoretical results and experimental 

data for Kjj are close to 1 up to 6c.ni. — 90°, but data 
sets at different energies seem to be inconsistent as they 
show different angular dependence. In contrast, theo- 
retical predictions at the three considered energies show 
nearly the same angular dependence for Ocm. < 90° . 

As already mentioned, above Ep — 20 MeV the min- 
imum of the elastic differential cross section is under- 
predicted. In order to establish the importance of the 
3NF as a means to cure this discrepancy we study the 
effect of the A-isobar excitation on both the differential 
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FIG. 4. (Color online) ^He analyzing power Aoy for elastic p-^He scattering at 7.03, 8.52j_^10.03, 13.6, 19.4, 21.3, 25.0, 30.0, 
and 35.0 MeV proton energy. Curves as in Fig. [T] The experimental data are from Refs. [30l. [32. l36l. [37l] . 
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FIG. 5. (Color online) p-'^He spin correlation coefficient Ayy 
for elastic p-^He scattering at 7.03, 8.52, 10.03, and 19.4 MeV 
proton energy. Curves as in Fig. [l] The experimental data 
are from Refs. [301.133]. 
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FIG. 6. (Color online) p-^He spin correlation coefficient 
Axx for elastic p-'^He scattering at 8.52 and 19.4 MeV proton 
ener gy. Curves as in Fig. [l] The experimental data are from 
Ref.ll. 
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FIG. 7. (Color online) Proton spin transfer coefficients K^ , K^ , and Ky for elastic p-^ He scattering at 
16.23 MeV proton energy. Curves as in Fig. [T] The experimental data are from Refs. [381. 1391]. 



.52, 10.77, and 



cross section and proton analyzing power at 30 MeV pro- 
ton energy. This has been done before at energies below 
three-cluster breakup threshold J4(|| and we follow here 
the same procedure. The results are shown in Fig. |S] in 
a way that one can single out the A-isobar effects of 2N 
nature, the so-called 2N dispersion, and of 3N and 4N 
nature, the 3NF and 4NF. The competition between 2N 
dispersion and 3NF, often found in the 3N system, is 
well seen also here for the differential cross section. As 
shown in Fig.|S]dispersive effect increases the discrepancy 
with data (dashed-double dotted curves) while 3NF and 
4NF effects reverse that trend for the differential cross 
section (dotted curves). Nevertheless, when the two ef- 
fects are put together the net result is an improvement 
towards the data (solid curves) but not quite enough to 
bridge the original gap. For Ay only the dispersive effect 
around the minimum is visible; it moves the predictions 
away from data. 



IV. SUMMARY 



In summary, we performed fully converged proton- 
■^He elastic scattering calculations with realistic poten- 
tials above the three- and four-cluster breakup threshold. 
The symmetrized Alt, Grassberger, and Sandhas four- 
particle equations were solved in the momentum-space 
framework. We used the complex energy method whose 
accuracy and the efficiency is greatly improved by the nu- 
merical integration technique with special weights. The 
pp Coulomb interaction was included rigorously using the 
method of screening and renormalization. 

The differential cross section exhibits rapid energy de- 
pendence and, in the minimum region around Oam. = 
120°, also sensitivity to the NN interaction model. The 
calculations using the INOY04 potential describe the ex- 
perimental data well up to Ep — 20 MeV but under- 
predict the differential cross section in the minimum at 
higher energies; other potential models fail even more. 
In contrast, most of the calculated spin observables show 
little sensitivity to the interaction model, and also the 
dependence on the beam energy is weaker than below 
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FIG. 8. (Color online) Differential cross section and proton 
analyzing power Ay for elastic p-"^He scattering at 30 MeV 
proton energy. Results obtained with CD Bonn (dashed- 
dotted curves) and CD Bonn + A (solid curves) potentials 
are compared with the experimental data from Refs. 26, 33]. 
In addition are shown the predictions including A effects of 
2N nature (dashed-double-dotted curves) and of 3N and 4N 
nature (dotted curves). 



the three-cluster breakup threshold. The overall agree- 
ment with the experimental data for the spin observables 
is quite good, considerably better than in the low-energy 
p-'^He scattering which is affected by P-wave resonances. 
In particular, the peak of the proton analyzing power 
Ay that is strongly underpredicted at low energies, is re- 
produced fairly well above E^ — 20 MeV but there is 
discrepancy in the minimum. The observed sensitivity 
to the NN interaction model seems to be mostly due to 
different predictions of the ^He binding energy; the cal- 
culations using the INOY04 potential with correct '^He 
binding provide the best description of the experimental 
data. 



We also studied the effect of three- and four-nucleon 
forces through the explicit inclusion of the A-isobar exci- 
tation. We found that A-generated many-nucleon forces 
significantly improve the description of the differential 
cross section but have almost no effect for the proton an- 
alyzing power. However, there are also quite strong dis- 
persive A-isobar effects that often reduce or even reverse 
the effect of 3NF. Therefore the total A-isobar effect, al- 
though beneficial, is not large enough to bridge the gap 
between the differential cross section data and calcula- 
tions. It even increases the discrepancy in the minimum 
of Ay . It might be possible that using the standard ap- 
proach of including static 3NF one might be able to ex- 
plain the data at higher energies, particularly using Effec- 
tive Field Theory generated interactions [4l|, |42| • Exten- 
sion of the method to other reactions in the four-nucleon 
system is in progress. 
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